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The properties of long gamma-ray bursts in massive 
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ABSTRACT 

Wc consider a popular model for long-duration gamma-ray bursts, in which the 
progenitor star, a stripped helium core, is spun up by tidal interactions with a black- 
hole companion in a compact binary. We perform population synthesis calculations 
to produce a representative sample of such binaries, and model the effect that the 
companion has on material that falls back on to the newly- formed black hole. Taking 
the results of hydrodynamic models of black-hole formation by fallback as our starting 
point, we show that the companion has two principal effects on the fallback process. 
First, a break forms in the accretion curve at around 10"* s. Secondly, subsequent to 
the break, we expect to see a flare of total energy around 10^° erg. We show that 
the break and flare times are set largely by the semi-major axis of the binary at the 
time of explosion, and that this correlates negatively with the flare energy. Although 
comparison with observations is non-trivial, we show that our predicted break times 
are comparable to those found in the X-ray light curves of canonical long-duration 
gamma-ray bursts. Similarly, the flare properties that we predict are consistent with 
the late-time flares observed in a sub-sample of bursts. 

Key words: Stars: gamma-ray burst: general; stars: binaries: close; stars:evolution; 
stars: supernova: general 
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1 INTRODUCTION 

The launch of the Swift satellite in 2004 has provided a 
massive step forwards in our abili ty to observe gamm a-ray 
bursts, in particular at early times (|Gehrels et al. 1120041 ). The 
combination of rapid triggering and localisation with pan- 
chromatic follow-up has provided a large and growing sam- 
ple of burst light-curves and spectra with which to confront 
theoretical models. Using this data, theoretical work con- 
tinues to elucidate the progenitors of gamma-ray bursts. In 
this paper we focus on long gamma -ray bursts (LGRBs), for 
which the model of lWooslevI (|l993l ) has stood up well to the 
last two decades of observations. This posits the origin of 
LGRBs as accretion of a stellar-mass quantity of material 
from an angular-momentum-supported disc around a newly- 
formed stellar-mass black hole. This model is supported by 
strong observational evidence for an association bet ween 
LGRBs and type Ib/c supernovae dStanek et al.l 120031): for 
a recent review of the evidence see lHiorth fc Bloom 1 201lf ). 
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The specific rates of gamma-ray bursts, even once cor- 
rected for a narrow beaming angle, are consistent with 
them forming from a small fraction of type Ib/c supernovae 
(,Podsiadlowski ct al..,2004i ). 



The main theoretical problem with this scenario is that 
a rapid rate of stellar rotation is needed at the moment of 
core collapse, in order to endow the core with sufficient an- 
gular momentum that its outer parts are unable to collapse 
and instead form the accretion disc. However, the presence 
of strong stellar winds in the Wolf-Rayet stars that lead to 
type Ib/c supernovae would be expected to brake the rota- 
tion of the stars and spin them down. A number of possi- 
bilities have been considered to overcome this problem. For 
single stars, the most promising option is low metallicity, 
which then leads to chemi cally homogeneous ev olution of 
rapidly-rotating stars (e.g. lYoon fc Langeill2005l ). The low 
metallicity suppresses stellar mass loss and thus prevents 
wind-driven spin-down. However, more recent observational 
evidence suggests that a significant fraction of LGRBs form 
in higher-metallicity environments where this mechanism 
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is no t plausible (e.g. iLevesgue et al.l I2OI0I : ISvensson et al.l 
|2012| ). 

Models proposed for the formation of gamma-ray bursts 
at higher metallicities mostly focus on exploiting the reser- 
voir of angular momentum stored in the orbit of a binary 
stellar system in order to provide the spin of the explod- 
ing star. Here we focus on angular momen t um transfer 
by tidal interaction JPodsiadlowski et al.ll2004l : llzzard et al.l 
l2004l : lLevan et al.ll2006l ). In this scenario, tidal locking of the 
progenitor core with a massive black-hole companion in the 
latest stages of its evolution spins the core up to rapid ro- 
tation rates. In Section[2]we use the BSE stellar population 
synthesis code to demonstrate the viability of this pathway 
and measure its frequency. We then analyse the effects of the 
binary companion on the accretion rate in the subsequent 
LGRB, in an attempt to confront this model with the ob- 
servational data. In Section[3]we describe our model, and in 
Section m our results. Section [5] contains a comparison with 
observed long-duration gamma-ray bursts, and Section |6] a 
summary. 



2 BINARY POPULATION 

To obtain a representative population of systems we syn- 
thesise a population of BH-BH binaries , utilising the rapid 
binary population synth esis code BSE (iHurlev et ani2002l ') 
with the modifications o f IChurch et al.l (|201ll ) which for con- 
venien ce we summarise here. We have modified BSE fol - 
lowing iBelczvnski et al] (|2002l ) and IBelczvnski et all (|2008l ') 
to include a more realistic prescription for compact ob- 
ject masses, hypercritical accretion during common envelope 
evolution, and delayed dynamical instability in mass trans- 
fer from helium stars in binaries with a lar ge mass ratio. 
We di stribute stellar masses according to the iKroupa et ahl 
(| 19931 1 IMF. Binary semi-major ax;es are chosen from a dis- 
tribution fiat in log a between 1 and 10* Rq . The initial ec- 
centricities of the binaries have little effect on their evolution 
and are set equal to 0.1 for all binaries. Neutron stars re- 
ceive a natal kick chosen from the bimodal distribution of 
lArzoumanian et al.l (|2002l ). We do not apply a kick to black 
holes upon formation. 

2.1 Criterion to obtain a tidal-spinup gamma- ray 
burst 

Following iLevan et all (|2006l '). we require that the material 
at the edge of the stellar core must have sufficient angular 
momentum to fall back into a disc at a radius outside the 
last stable orbit. This means that the specific angular mo- 
mentum j of material at the edge of a core of mass Mc must 
satisfy j > \/&GMc/c. Taking the core to be rotating as 
a solid body and be tidally locked to the binary, we apply 
Kepler's law to obtain a critical orbital separation Ocrit of 

where we have taken a typical core radius for massive stars 
of -Rc = 0.2 J?0 and Mt is the total mass of the binary at 
this point. For binaries that satisfy this criterion we assume 
that tidal locking takes place at the point where the orbital 
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Figure 1. The population of binaries that we expect to produce 
tidal spin-up gamma-ray bursts. In the bottom-left panel we plot 
the mass of the companion black hole, Mi , against the binary sep- 
aration before the second supernova, a. The top-left panel shows 
the mass lost in the supernova, AM2, against Mi. The lower-right 
panel shows a against AM2. In BSE AM2 is a function of M2 
only. Red dots show all the binaries produced by our population 
synthesis. Blue squares show the subset of 17 representative sys- 
tems that we chose to investigate and the filled blue square shows 
the location of our typical system. 

separation is minimised. This gives the specific angular mo- 
mentum of the material at the edge of the core as 

A„,=6,97X10'>.„.^-^(|=|)-"', ,2) 

For cores that match this criterion, the material further 
inside the core also acquires a substantial angular momen- 
tum. The black hole that forms from such a core will have a 
dimensionless spin parameter close to unity. This makes it 
a good candidate for powering a long gamma -ray burst via 
the mechanism of iBlandford fc Znaiekl (|l97iD . 

2.2 Derived population 

The population of binaries that we obtain from our BSE 
runs is shown in Figure [T] We plot the relationships between 
the companion mass. Mi, the pre-supernova binary semi- 
major axis, a, and the mass lost by the exploding star, AM2. 
The binaries where the criteria for a tidal-spinup LGRB are 
satisfied are marked with red dots. We choose a sub-sample 
of 17 systems that cover the parameter space evenly to in- 
vestigate further; these are marked with blue squares. 

To calculate a repr esentative rate for the s e syst ems we 
adopt the approach of IPodsiadlowski et all (|2004l '). along 
with their rates for LGRBs and core-collapse supernovae 
in a typical galaxy. We present the rates for our binaries in 
Table [1] To calculate our rates we employ a iKroupa et al.l 
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Object Rate 

yr-l/galaxy 

Core-collapse supernovae 7 X 10"'^ 

Gamma-ray bursts (beaming angle of 5°) 3 X 10" '' 

Tidal spin-up LGRBs (conservative estimate) 2.6 X 10~^ 
Tidal spin-up LGRBs (optimistic estimate) 1.2 X 10~^ 



and subsequently falls back on to the black hole receives the 
same kick as the newly-formed black hole itself. This is the 
expected behaviour if the kick is caused by a mechanism in- 
trinsic to the supernova explosion, which previous work has 
found to be the mos t likely case for the G alactic black-hole 
binary Cygnus X-1 (|Axelsson et al.|[201ll ). 



Table 1. Estimated event rates for a typical galaxy. Values for 
core-collapse supernovae a nd gamma-ray bursts are taken from 
iPodsiadlowski et al.l | |2004) . 



IMF and a minimum mass for a core-collapse super- 
nova of 8 Mq . Our conservative estimate takes a binary frac- 
tion of 0.5 for all stars and assumes that their masses are 
drawn independently from the IMF. The more optimistic 
estimate - which in our judgement is more realistic - takes 
the binary fraction to be unity for massive stars and em- 
ploys a distribution of mass ratios f{q) = 2q, which biases 
the sample towards equal-mass binaries. 

Our rates, taken at face value, are quite small; even the 
optimistic estimate is lower than the LGRB rate. This is 
relatively unsurprising; we require a binary that forms two 
black holes, one of them rather massive, in a binary with an 
unusual configuration. It is perhaps unsurprising that such 
binaries are relatively rare compared to, for example, NS- 
NS binaries. However, as we show, making a simple change 
to the mass-ratio distribution in the binary can make a very 
large difference to the predicted rate. We have not explored 
further any of the uncertainties in our assumptions about the 
binary evolution and how they might affect the population of 
binaries produced, but we are confident that, were we to do 
so, we would be able to match the observed GRB rate using 
reasonable assumptions. For example, the number of suffi- 
ciently spun-up binaries that we produce is rather sensitive 
to our assumptions about the radius of the helium core at 
the closest point in the evolution. We have also not evaluated 
whether the injection of angular momentum into the falling- 
back material by interaction with a companion can increase 
the population of binaries with sufficient angular momen- 
tum to form a disc. There are also potentially other possible 
routes to making gamma-ray bursts via tidal spin-up that 
we ha ve not considered. For example. iMoreno Mendez et al.l 
(|201ll ) show that the high spin of black holes in two Galactic 
binaries is consistent with t hem being the r elics of gamma- 
ray bursts, as predicted bv lLee et al.l (|2002l ). Such binaries 
are likely to be relatively common compared to the massive 
double black hole binaries that we consider here. 

2.3 The effects of natal kicks 

Ther e is some evidence that black holes receive n atal kicks 
(e.g. iBrandt et all 1 19951 : iGualandris et all |2005| . see also 
Repetto et al. submitted). To investigate the likely effect of 
natal kicks on the processes that we are modelling we con- 
sider the effect of a moderate kick - 100 km s^^ - in the four 
directions in the orbital plane: ±{x,y}. We do not consider 
kicks in the z direction because this breaks the symmetry of 
the problem about the orbital plane, and hence evaluating 
the behaviour of the subsequent accretion flow is more com- 
plex. In each case we assume that the material that is ejected 
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3 MODELLING THE FALLBACK PROCESS 

To treat the fallback pr ocess we utilise the simulations of 
iMacFadven et all (|200ll ). They find a broadly self-similar 
behaviour of the accretion of material on to the nascent 
black hole , with an early p lateau phase followed by the t~^^^ 
decline of IChevalierl (| 19891 ). We fit this with a two-piece so- 
lution: 

100 S < f < fplat eau /n\ 

-5/3 i-^i ^ ' 

^ ^ '-plateau 

where logi(,(tpiateau - 100 s) = 2.2; i.e. fpiatoau 260 s, and k 
is chosen to obtain the correct total fallback mass. To obtain 
the fallback mass we assume that the inner 2 Mq of the final 
black hole mass is acquired during the initial core collapse; 
the remainder forms the fallback mass. This prescription 
leads to a fallback mass of between 1 and 3 Mq . 

3.1 Particle trajectories 

By using the results of hydrodynamic simulations of black- 
hole formation in fallback supernovae in the form of Equa- 
tion[3]we avoid having to compute the hydrodynamics of the 
accretion shocks ourselves. Instead we treat our particles as 
interacting only under gravity. We launch our particles with 
a distribution of velocities, chosen so that, in the absence of 
a companion star, the particles would fall back with times 
that reproduce the accretion rate given by Equation |3] For 
purely radial motion the return time troturn for a particle 
launched at radius ro from a star of mass Mi is given by 



^return — 



GMi Vl-Q 



a _i / a 
■ tan 



1 - a 



1 - Q 



(4) 

where a is related to the particle's launch velocity vq ac- 
cording to 



2GMi 



ro 



(5) 



We invert Equation. [3] numerically in order to produce a 
distribution of a from the accretion rate of Equation. [S] 
Because the star is rotating the initial velocities are not ra- 
dial; however, because ro is small compared to the maximum 
distance that the particle reaches from the exploding star. 
Equation |4] suffices. We choose the value of ro so that the 
rotational velocity at the equator is 20% of the escape veloc- 
ity at that radius; thus even at launch the particle velocities 
are predominantly radial. 

We follow the motion of the particles under the grav- 
itational force of the two black holes. As we are using 
Equation |3] to account for the self-interaction of the gas 
we ignore the mass of the gas particles and any non- 
gravitational forces; in practice, therefore, we solve the re- 
duced three-body problem numerically. We use the Bulirsch- 
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Stoer method as given bv lPress etHI j 19921 ) to integrate the 
orbits of the stars and particles. 

We follow the motion of the particles until they cross 
the orbital plane {z = 0) for the first time. The angular 
momentum vector of the stellar rotation is parallel to that 
of the orbit as it is induced by tidal spin-up of the explod- 
ing star: hence the system is symmetric under reflection in 
the orbital plane. Because we can expect any given particle 
that crosses the orbital plane to have a counterpart with the 
equal and opposite z-velocity a disc will form in the plane. 
We measure the mass and momentum flows into the disc 
by monitoring the rate at which particles cross the plane 
and the specific angular momentum that they possess when 
they do. We assume that any particles that cross the plane 
sufficiently close to a black hole to be within the last stable 
orbit are accreted directly. 

3.2 Accretion disc 



Following IPerna et al] (|2006h we use a simple treatment of 
the accretion disc based on the viscous timescale. We assume 
that a ring of material, represented by one of our particles 
that has accreted into the x — y plane, will move inward 
through the disc via viscous dissipation of angular momen- 
tum. The rate at which a particle moves inward is given 
by 

dr _ r 
di" ^ t^' 

with characteristic viscous timescale to given by 



(6) 



to = 



^^dis 



Kop 



(7) 



for a disc of scale height H and radius 7?. The orbital angular 
frequency at radius r in the disc around the compact object 
with mass Mi is given by 



Kcp 



GMi 



(8) 



We take the ratio -Rdisc/^^disc = 10 to be fixed through- 
out the evolution. Although the disc is probably less fiat 
than this at early times the viscous timescale then is short 
enough for the results to be basically unaffected. A particle 
unaffected by other particles, arriving at the disc at time r, 
then spirals towards the central black hole according to 



r[t) = r(r) 



3t 



-1 2/3 



2fo(r) 



(9) 



We set r(r) to be the particle's circularisation radius; that 
is, we assume that the material that it represents will self- 
collide to form a circular ring, conserving specific angular 
momentum. Thus the initial radius at which the particle 
joins the disc is given, in terms of its specific angular mo- 
mentum j, by 



r(r) = r-circ = 



(10) 



In the process of circularising, i.e. moving between the radius 
from the central object at which it enters the plane and that 
at which it circularises, the material may collide with gas 
already present in the disc. If this happens we merge the 
two colliding particles and calculate a new specific angular 



momentum, and hence circularisation radius. This process 
is repeated until the particle has successfully circularised. 
Particles that circularise outside the instantaneous Roche 
lobe of either star are assumed to be lost from the system. 
We have tested that the accretion history resulting from this 
treatment converges as we increase the number of particles, 
and hence the mass resolution of the simulation. 



4 RESULTS 

We present first the results for a typical binary, indicated by 
the filled square in Figure [T] It has a companion black-hole 
mass of Ml = 8.24 Mq; the newly- formed black hole has 
a total mass, including the material falling back, of M2 = 
4.29 Mq, having lost AM2 = 3.49 M© in the explosion. The 
semi-major axis of the orbit at the time of explosion is a = 
4.01 i?Q. 

In Figure [5] we plot the location of each particle in the 
simulation when it crosses the orbital plane. The exploding 
star is initially moving in the -\-y direction. As it proceeds 
round its orbit, the particles that fall back towards it have 
been affected to a greater degree by interaction with the 
companion, so their deviation from the stellar position is 
increased. These particles typically also have more specific 
angular momentum and hence enter the accretion disc at 
larger radii. After a some time - about one quarter of the 
post-explosion orbital period of the binary - the particles 
start to mostly fall back outside the Roche lobe of the star 
and hence the mass inflow rate into the disc reduces sharply. 
Some particles instead, however, fall into the Roche lobe 
of the companion black hole (blue dots in Figure [2]) . The 
accretion rate on to the two stars as a function of time in 
this typical model is shown in Figure [3] 

Our typical model shows a number of features that are 
common to the models that we have run. There is an ini- 
tial plateau followed by a declining power-law accretion rate. 
This follows directly from our input accretion rate; the par- 
ticles falling back at early times are largely unaffected by 
the presence of the companion. At roughly 10** s there is a 
break in the accretion curve, which transitions to a steeper 
power law. Finally there is a period of zero accretion at 
about 1.75 X 10'* s, followed by a flare that lasts for several 
lO'' s. During this flare there is additionally some accretion 
on to the companion black hole. 



4.1 Fits to the accretion behaviour 

To gain some insight into the implications of our model we 
fit each of our 17 characteristic systems with a broken power- 
low model of the form 



M oc 



t^l 200S < t < ibroak 

<t <to, 



(11) 



where 71, 72 and fbrcak are obtained by fitting. The time at 
which the accretion first ceases, to, is obtained from the ac- 
cretion curve. In Figure 3] we present the results of the fits. 
The bottom-left panel shows the break time, tbreak, plotted 
as a function of 71, the power-law index of the portion of 
the accretion curve before the break. The plot shows that 
the accretion curves are somewhat steeper than the t^^^'^ in- 
put power-law; shorter break times are correlated with more 
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Figure 2. The locations where particles cross the orbital plane 
in our reference simulation. The large red and blue points are the 
initial positions of the exploding star and companion black hole 
respectively. The red and blue solid lines show the post-supernova 
orbits of the two stars. Initially the exploding (red) star is moving 
in the +y direction and the companion black hole (blue) in the —y 
direction. The red dots are particles that fall into a disc around 
the newly-formed black hole; similarly, the blue dots are particles 
that fall into a disc around the companion black hole. Grey dots 
are particles that fall further from either object than its Roche 
lobe radius and hence are expected not to be accreted. The black 
dashed lines show the initial radii of the Roche lobes around the 
two stars. The origin of the co-ordinate system is the centre of 
mass of the post-explosion system. For clarity only 1% of the 
particles used in the final simulation are plotted. 



steeply declining accretion. The power-law index of the sec- 
ond, post-break power law segment is much steeper, typi- 
cally between -5 and -10, and does not show any correlation 
with the other parameters. 

In Figure[5]we present the correlation between the semi- 
major axis of the binary, a, and the accretion break time. 
The evident correlation has a simple explanation; a wider 
orbit means that the material has to travel further before 
it interacts significantly with the companion black hole, in- 
creasing f break- This in turn also means that the angular mo- 
mentum imparted to the material which is accreted before 
the break is lessened; hence there is a positive correlation 
between a and 71. 

In both Figures|4]and[5]we plot both systems computed 
under the assumption of black-hole kicks, as well as unkicked 
systems. We find that the behaviour of the kicked and un- 
kicked systems is broadly similar to the unkicked case, and 
there are no marked trends between the different kicks. This 
suggests that a moderate kick is unlikely to have a strong 
effect on the fallback process. 

© 2012 RAS, MNRAS 000, [iHS] 
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Figure 3. The accretion rate on to the black holes in our typical 
example system. The grey dashed line shows the accretion rate 
of Equation [S] which would be obtained in the absence of a com- 
panion. The red solid line shows accretion on to the newly-formed 
black hole (star 2) whilst the blue dotted line is for the companion 
(star 1). The plot includes the effect of the accretion disc and the 
Roche lobes of the two stars. Of the 2.29 of material in the 
fallback shell 1.90 is accreted on to star 2 and 0.02 on to 
star 1; the remaining 0.37 is lost from the system. Character- 
istic features are (a) the initial plateau and roughly t~^/^ decline 
phases, (b) the break in the accretion rate at roughly 10* s fol- 
lowed by a steep decline, and (c) some late-time accretion activity 
(at around 3 X lO'^s). The right-hand axis shows the energy in- 
jection from accretion, assuming that the efficiency of conversion 
of mass to energy is »y = 0.1. 



5 COMPARISON WITH OBSERVATIONS 

It is difficult to compare our results directly with observa- 
tions. In order to go from the black hole's accretion history 
to the observed photometry and spectral features it is neces- 
sary to model the propagation of the fireball and jet through 
the exploding star, a highly involved procedure. What we 
can say, however, is that we expect to see evidence of con- 
tinued energy injection during the first approximately 10* s 
in the rest frame of the accreting black hole. We predict 
that the rate of this energy injection would decline quite 
strongly, although we should caution against the use of our 
fitted values of 7 as the decline of the injected luminosity. 
As the system evolves and the accretion rate decreases the 
accretion morphology will become more disc-like and less 
advection dominated. This will increase the fraction of the 
accretion luminosity available to power a jet at later times. 

Secondly, we predict a sharp reduction in the accretion 
luminosity at times of around 10* s, which we would expect 
to appear as a break in the light cur ve. These times are 
consistent with typical times found by lEvans et al.l (|2009l ) 
for the end of the plateau phase in canonical light curves, 
whose distribution also peaks at about 10* s. Light curves 
that show a single, steepening break also have break times 
around 10* s, although the distribution is much wider. Our 
break times are also consiste nt with the sub set of the sample 
of light curves presented bv lLi et all (|2012l ) that show shal- 
low decay features. Interpretation of the optical light curves 
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Figure 4. Correlations between different parameters in our 
power-law fits to tiie accretion rates on to tlie two black lioles. Tfie 
power-law applicable to the decay between 200 s and the break 
time tbreak is denoted by 71, whereas the power law between 
*break s-^d the first cessation of accretion is denoted by 72- The 
bottom- left panel shows the correlation between 71 and tbreak i 
the top-left panel that between 71 and 72, and the bottom-right 
panel that between 72 and fbrcak- Black crosses are for models 
without natal kicks, brown filled circles and squares are models 
with kicks in the +y and —y directions, and green open circles 
and squares are models with kicks in the -l-a; and —x directions. 
Only the first of these three shows a significant correlation. 
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Figure 5. The correlation between the binary semi-major axis 
a and the break time tbreak- Wider binaries only affect particles 
falling back at later times, thus making tbreak later for larger 
semi-major axis. Black crosses are for models without natal kicks, 
brown filled circles and squares are models with kicks in the +y 
and —y directions, and green open circles and squares are models 
with kicks in the +x and —x directions. 



is complicated, however, by spectral breaks which can occur 
at relevantly late times. 

The physical origin of the break in our models is the 
truncation of accretion by the Roche lobe of the exploding 
star. Material that falls outside the Roche lobe is not bound 
to either star, and hence will be lost from the system. As 
time progresses, the material falling back has traversed a 
larger distance from the exploding star and hence is more 
affected by the gravitational force of the companion. Thus 
it falls back at a greater distance from the exploding star. 
The break time can be understood as the time at which 
this distance is the Roche lobe radius of the star. This is 
the origin of the correlation visible in Figure [5] as in wider 
binaries the Roche lobe is larger and thus the break time is 
later. 



5.1 Late-time highly-variable accretion 

The majority of our models include some form of late- 
time highly-variable accretion. This typically occurs at a 
few times tbreak, after the second power- law component has 
ended. In our typical system, for example, after the main 
accretion episode has ended at around 2 x 10* s there is a 
renewed episode of accretion peaking at around 3 x 10* s. 
Such activity might power late-time flares. In Figure [S] we 
show the properties of the flares in those of our models 
that show them. We assume that, once the accretion en- 
ergy conversion efficiency has been taking into account, all 
the energy liberated goes into the flare. Having done this, 
the luminosity generated is very comparable to that seen in 
late-time flares observed in the X-ray and optical. For com- 
parison, we plot the energy in the late-time flares of three ob- 
served long gamma-ray bursts, GRB 050502B, GRB 070107 
and GRB 070318. This i s the subset of the sample con- 
sidered by ICurran et al.l (|2008l l that have reasonably con- 
strained photometric or spectroscopic redshifts and which 
are l ong bursts. We ta ke the light curve s for these flares 
from [ Evans et all (|2009l ) and redshifts from lXiao fc Schaeferl 
(|20lj r Although the observed flares appear to be somewhat 
more energetic than our predictions, it is only necessary for 
them to be beamed to one tenth of the sky at the time 
of the flare to be consistent with our predicted values. We 
re-produce the observed decline i n flare energy with time, 
which is also seen in optical data (|Li et al.ll2012l V 



6 DISCUSSION 

We have considered tidal spin-up compact binaries as a pos- 
sible source for long-duration gamma-ray bursts. We show 
that, with appropriate assumptions, a population synthesis 
approach predicts that these binaries will be produced. Fur- 
thermore we show that, given two different - but plausible 
- assumptions about the distribution of mass-ratios of high- 
mass binaries we predict two very different formation rates. 
Hence we suggest that the power of population synthesis 
calculations to constrain progenitor model s by analysi s of 
expected rates is rather weak, as argued bv lToutI |200i) for 
the similar problem of Type la supernova progenitors. 

A potential means of distinguishing this progenitor sce- 
nario, however, is that in the case of a tidally spun-up pro- 
genitor, the star exploding as a type Ib/c supernova has 
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Figure 6. Episodes of highly-variable late-time accretion in our 
models. Wc plot the total mass accreted in late-time flares, Mfiarci 
as a function of the break time tbroak- As in previous figures, black 
points are for models without natal kicks, brown filled circles and 
squares are models with kicks in the +y and —y directions, and 
green open circles squares are are models with kicks in the +x 
and —X directions. Flare energies (right-hand axis) are calculated 
assuming mass is converted into energy with an efficiency of = 
0.1. The pink upper limits represent the energy in late-time flares 
from three long-duration gamma-ray bursts, (a) GRB 050502B, 
(b) GRB 070107 and (c) GRB 070318. In calculating these points 
we assume that the energy emission at the time of the flare is 
isotropic. 



a black-hole companion. We have modelled the effect of 
the companion on the fallback process to investigate what 
observable consequences it would have for the gamma-ray 
burst. We conclude that our models possess the following 
general features: 

• Before about 10* s, an accretion rate onto the newly- 
formed black hole that declines as ^-i ssio.os^ 

• a break in the accretion rate at around 10* s, with a 
steeply declining accretion rate thereafter, and 

• in the majority of cases, highly-variable late-time ac- 
cretion, leading to flares at times between 10* and 10^ s, of 
total energies around 10^° erg. 

These general features are broadly consistent with the 
observations of long-duration gamma-ray bursts. Bursts 
with a canonical-type light curve show evidence of energy 
injection up until roughly the times that we predict, with 
a break there a fter. F or the three bursts from the sample of 
ICurran et al.l ()2008h that have a well-deflned redshift and 
late-time bursts, the properties of the bursts agree with our 
predictions, assuming a reasonable beaming fraction. We re- 
produce the tentative correlation between the time of the 
flare and the energy. We caution, however, that our model 
on its own cannot explain flares that occur prior to 10* s. 
It does not preclude such flares from occurring; however, 
they must have a separate origin, such as instabilities in the 
accretion flow onto the newly-formed black hole. 

Although we have only considered tidal spin-up by a 
massive compact companion, it is also possible that the ex- 
ploding star could be spun up by a main-sequence compan- 



ion. iLee et al.l (|2002l ) consider the evolution of the progeni- 
tors of low-mass X-ray binaries. They show that tidal forces 
from a low-mass main-sequence companion and frictional 
effects during common-envelope evolution can be sufflcient 
to spin up the outer layers of the helium star, which subse- 
quently collapses to give a black hole. We expect gamma-ray 
bursts produced by such a system to also show breaks and 
flares similar to those described in this paper; the binary has 
a similar semi-major axis and the same processes will oper- 
ate. We intend to undertake a quantitative study of these 
binaries in a subsequent work. 

Finally, we investigate the effect of black-hole kicks on 
our scenario and show that a moderate kick is expected to 
have a relatively small effect on the observed properties of 
the burst. We caution, however, that our model can only 
encompass kicks that are directed within the orbital plane 
of the binary. Kicks with a component parallel to the an- 
gular momentum vector would disturb the geometry of the 
situation; this could make the accretion flow more disc-like 
and less spherical, which we would expect to have significant 
consequences for the baryon loading of the resulting burst. 
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